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The Abelian Higgs model with or without external particles is considered in curved 
space. Using the dual transformation, we rewrite the model in terms of dual gauge fields 
and derive the Bogomol'nyi-type bound. We find all possible cylindrically symmetric vortex 
solutions and vortex-particle composites by examining the Einstein equations and the first- 
order Bogomol'nyi equations. The underlying spatial manifold of these objects comprises a 
cylinder asymptotically and a two sphere in addition to the well-known cone. 
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I. Introduction 

A classic example of solitons in field theories is the topologically stable vortex solution 
in Abelian Higgs model [p], J2[. Such vortex-like excitations have been used to describe the 
system of flux-tubes in type-II superconductors JT|. These objects have played consider- 
able roles also as cosmic strings. For example, phase transitions in early universe give rise 
to string-like defects and they generate the density fluctuation which leads to a possible 
explanation of galaxy formation [H |J. 

In this paper, we consider U(l) local strings of Einstein Maxwell Higgs theory. However, 
the model may be too complex to handle analytically if one considers arbitrary shaped cosmic 
strings coupled to gravity. Thus one may adopt local cosmic strings which are infinitely 
straight in one direction and remain in equilibrium, which describes an idealized but a 
physically relevant situation. This assumption simplifies the system to that of Nielsen-Olesen 
vortices coupled to Einstein gravity in (2+1) spacetime dimensions, and the zero size limit 
of such objects can be identified with point particles on a plane ||. Then (2+1) dimensional 
Einstein Maxwell Higgs theory with or without the coupling of external point particles is 
of our interest, and we examine possible vortex configurations whose centers point particles 
may stick to. Furthermore, if we choose a specific form of scalar potential, static local 
vortex configurations can be solutions of a first-order differential equations, which satisfies 
second-order Euler-Lagrange equation automatically |6|, 0]. These Bogomol'nyi equations 
have been established in the theory coupled to gravity despite the difficulty of constructing 
the gravitational energy p, [|, [[0| . For such vortex solutions to this Bogomol'nyi equation, we 
show that they are exact solutions of Einstein equation and that the cosmological constant is 
zero for a general stationary metric and arbitrary distribution of point particles. Specifically, 
the Bogomol'nyi- type bound is saturated when the metric is static ||. 

Another way of understanding the role of vortices in U(l) gauge theory is to recapitulate 
it in dual formulation and elicit the physically relevant aspects such as phase transitions 
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induced by vortices and classical dynamics of vortices [JTTJ, [12| . Here we derive the dual- 
transformed version of Einstein Maxwell Higgs theory in (2+1) dimensions, and show that 
the Bogomol'nyi-type bound can be attained within this formulation. It is found that the 
extended objects corresponding to the vortices in the original theory do not carry dual 
magnetic flux but dual electric field, and the energy which is proportional to Euler invariant 
is expressed by the sum of the spatial integral of electrostatic potential and the total mass 
of point particles in the Bogomol'nyi limit. 

Bogomol'nyi equations being set up, we investigate cylindrically symmetric configurations 
and find all possible vortex-particle solutions by explicitly proving their existence. Interest- 
ingly, there are solutions whose two manifold constitutes a cylinder in its asymptotic form 
or a two sphere, in addition to those solutions whose two manifold constitutes a cone. 

The rest of this paper is organized as following. In section 2 we introduce the model 
and derive Bogomol'nyi-type bound in a more detailed way by solving Einstein equations 
under general stationary metric. Section 3 is devoted to the dual transformation of Einstein 
Maxwell Higgs theory within path-integral formalism and to the study of various aspects 
of the dual formulation, including the derivation of Bogomol'nyi-type bound. In section 4 
we analyze cylindrically symmetric solutions and their global geometrical structures. We 
conclude in section 5 with some remarks. 

II. Einstein Maxwell Higgs Theory and Bogomol'nyi Bound 

We consider the system composed of vortices in Abelian Higgs model and the massive 
point particles in the presence of gravity. Einstein Maxwell Higgs theory coupled to a set of 
N external particles in (2+1) spacetime dimensions is described by the action 

*J ^gravity ~T iJmatter "T ^point particle 

= Jd 3 x^l [ - Y ^(R + 2A)- 1 -g^g^F tlu F pa + ±g^D^D v (f> - V(|0|) } 



_N_ r oo / d,T a rt r r a 



dx?, dx?, 
m a I ds\j g^ 1 

a=\ 

where = e* n |0|, D^<p = (d^ — ieA^cp, and A the cosmological constant. 
Equations of motion read 

-± = d v ( y /gF>»') = ef (2.2) 

±DAWM) = ~% (2-3) 

fl^-j9^(fl + 2A) = 8jrGT„„, (2.4) 



where j M is the conserved U(l) current, 



f = -\<r(j>D v <j>-D v <M>), (2-5) 



and Tn U is the energy- momentum tensor, 



1 A Z" 00 , „, , , ,dx"dx?, 



-"matter 



+— £ ™<* / ds a 5(x a (s a ) - x )—a—a.. (2.6) 

yfi 1 a=l °° a a 



We are interested in static soliton solutions of the equations of motion, specifically the 
Nielsen-Olesen vortices which are electrically-neutral and are characterized by the magnetic 
flux 

<$> = - J d 2 x l -^F l3) (2.7) 

where e*- 7 two dimensional Levi-Civita tensor density of e 12 = e 12 = 1. With an appropriate 
choice of the potential V(ip), let us look for multivortex-particle configurations in curved 
spacetime which satisfy Bogomol'nyi-type equations. If we assume that the vortices satisfy 
such equations, Einstein equations will turn out to be solved under the general stationary 
metric of the form 

ds 2 = N 2 (dt + KM) 2 - lijdx'cy, (2.8) 



where the functions N(x), Ki(x), and jij(x), (i,j = 1, 2), are independent of time. 

Under A = gauge condition, Oz-components of energy-momentum tensor T* vanish, 
so the corresponding solution of Einstein equations is 

K„ = *$=, (2.9) 

where K^ = diKj — djKi, and k is a undetermined constant. Taking N — 1 and rearranging 
the terms in the spatial integration of the 00-component of energy-momentum tensor, we 
have 

/ dfay/g Too 



+-7 y (A0 =F zV7Qfc7 W A0)(^0 =F »V7^7 m ^n^) 
±^.(V2T-|(i0| 2 -t; 2 ))| (2.10) 



2 N 

Z o=l 



Now we choose the scalar potential as 

^ = ^(|0| 2 -^ 2 ) 2 , (2-11) 



then the first two terms in the braces of Eq. ( |2.1C| ) are nonnegative definite and we obtain 
the so called Bogomol'nyi bound. The bound is attained for configurations satisfying the 
Bogomol'nyi equations 

A0TV7%7 ife £W> = O (2.12) 
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-F y = ±V2V. (2.13) 



In this Bogomol'nyi limit, the ij-components of energy-momentum tensor vanish and then 
Einstein equations force both an integration constant k and the cosmological constant A to 
be zero. Thus we have no off-diagonal components of metric Ki up to a gauge. The 00- 
component of Einstein equations relates the topological quantities, i. e. the geometric part 
gives the Euler number and the matter part, the sum of magnetic flux from vortices and 
total mass from point particles 

1 r v 2 N 

d 2 xJ^ 2 R=—\e§\ + Ym a . (2.14) 

16nG J V ' 2 ' ' ^ K ' 

After eliminating the gauge field A4 by use of a Bogomol'nyi equation ( |2.12j ) 

A t = ~($ft T Vl^l 3k dk In |0|), (2.15) 

and fixing the gauge for the spatial components of metric 

lij = -5 ij b(x i ), (2.16) 

we solve 00-component of Einstein equations 



''(•'■' > = < "' i]+m f TTn ,, iZ^, ,m. ) . ( 2 - 17 ) 



p e -(f-l) 

I [" I •■ _ z \ 2 U N \z — z I 2 ™" 

,llp=l l- 6 *p\ lla=l I* z a 



where h(z) (h(z p )) is a holomorphic (antiholomorphic) function and the tilded variables are 
dimensionless quantities 



,1 :.2 _ . .., J ;..2, , ... M n _ a -^..2 ^ ... ™>° _ 



z = x l + %x l = evfx 1 + ix z ), f = izi 1 G = 4ttGV, m a = — ^. (2.18) 

v nv 2 



If the spacetime of the vacuum (no particle (N = 0) and no vortex (n = 0)) is to be 
Minkowski, the harmonic function should be chosen to be zero, h(z) + h(z) = 0. Substituting 
Eq.( |2.17| ) and Eq. fl2.15| ) into Eq. ( |2.13| ), we obtain a single equation for the amplitude of Higgs 



field 



d 2 \nf 2 = e h+h [— — J - N , — (/ 2 -l)+47rV^(f-5 p ), (2.19) 

\llp=l l z n>l llo=l K z «l / p=l 



where d 2 is a Laplacian in flat two-dimensional space. Returning to the expression of Euler 
number in Eq.( |2.14j ) and inserting Eq.( |2.17] ) into it, we have 
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d 2 xd 2 \n\(t)\ 2 + fd 2 xd 2 



(2.20) 



From the above expression one may notice that, in addition to the contribution from the 
first two terms for ti ^ or m a ^ 0, there can exist contribution from the second term 
if the Bogomol'nyi equation (|2.19| ) contains the finite energy solution which behaves as 
|0| ~ \x\~ £ (s > 0) for large |af|. We shall show that it is indeed the case and present the 
detailed analysis for the existence of such solutions in section 4. 

III. Dual Formulation 

In this section we reformulate the theory by use of the dual transformation and re- 
derive the Bogomol'nyi limit in this formulation. The path integral for the system we are 
considering is given by 



F\e mT \I 



[dg^ v ] [dA^ 
x expi< / d 3 x y/g 



_L_ (jR + 2A ) - lgWg™F,„F p<T + \gTdM\h 



+ 22H0| 2 (<9^ - eA,){d v n - eA v ) - V{\4>\) 



N ,00 / 

+ Y, m a ds\g^ 
0=1 J -°° v 



dx^ dx% 
ds ds 



(3.1) 



Introducing an auxiliary vector field C M , we rewrite the interaction term between scalar 
field and gauge field as the following 

exp {tjd 3 x^ i<r|0| 2 (c^ - eA^duil - eA,)} 

= J[dC»]ll^expifd 3 x^L^C»C„ + g^C^d„n-eA u )\. (3.2) 

Classically the auxiliary field C M is nothing but the conserved U(l) current. 

Since the theory contains vortex and antivortex configurations, the phase of the scalar 
field need not be single-valued and then can be split into two parts 

Q(t,x) =Q{t,x)+r]{t,x). (3.3) 

The first term which describes a configuration of vortices and antivortices is defined by a 

multi-valued function 

x — x (t) 

e(t,f) = E(T)tan- 1 ^ , (3.4) 

where (xUt), xi{t)) is the position of a (anti-)vortex and (=f) denotes —1 for vortex and 1 for 
anti-vortex. The single-valued function 77 represents the fluctuation around a given vortex 
sector. Hence the path integral measure is divided into two contributions 

[dSl] = [dQ}[d V ], (3.5) 

i. e. the one for the sum over single- valued fluctuation around a given configuration of vortices 
and the other for that over all possible configurations of vortices, including the annihilation 
and creation of vortex-antivortex pairs. 
After ^-integration, 

/[*7]exp{i j d 3 xy/g g^C^n} » -U(V M C"), (3.6) 

we can rewrite a part of the path integral by introducing dual gauge field H^, 

J[dC,]^-5(V,Cn ■■■ = J[dH,][dC,]5(^C^ - X -& vp d v H p ) ■■■ . (3.7) 



Substituting Eq. (|3.7| ) into the path integral in Eq. (|3.2| ) and integrating out the auxiliary 
field C p , we obtain an effective theory described by the action 

S' = Jd 3 x^L-^(R + 2A)-^F ll F u + ^d fl \<f ) \dM-V(\<f ) \) 

1 1 & vp 1 1 

jg^sTH^H^ + ~—H„ u (-d p G - A p ) (3.8) 



4e 2 |(/»r ^ P 2y/g ^e 



A r°° , / dx a u dxl 
ds ds 
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where H uv = dnH u — d v H u and F^ = F up , the dual of field strength tensor. For per- 

forming A p integration, let us regard A p and F^ as independent variables and rewrite the 



path integral measure by introducing another auxiliary field N f 

J[dA,][dF»]5{F» - —d u A p ) ■ ■ ■ 



V9 
= J[dA,}[dF»][dN^ W^fg- exp U J d 3 x^g N^F» - ^d u A p )\ ■■■ . (3.9) 

Putting the above equation ( |3.9|) into the path integral and then integrating F^ and A^ 
fields out, we have 

J[dA p ] exp i J d 3 x ^g^-lg^p^ + H p F^ 

= J[dN p ] \{g^ 8{e^d v N p ) expzj d\jg 9 -^-{H p + N p )(H u + N„) . (3.10) 

The delta functional in Eq.( p.lO|) implies N p = —d P x f° r a single- valued scalar field x- 
Finally the path integral of dual transformed theory becomes 

Z = J\gdg^][dH^[\<f>\-^<f>\][de][d X ]exp{iS D }, (3.11) 

where the action of dual transformed theory is 

Sd = Jd 3 x^l [ - T ^(R + 2A) + ^d fl \<j ) \d u \ ( f ) \-V(\<j ) \) 



1 1 & vp 1 

g^g™H, v H pa + -—H^B + -g^(H, - d pX )(H u - d uX ) 



4e 2|0| 2 - - ~r~v» 2ey/g ^ " 2 

+ S^/-1 , '4*' (3 ' 12) 

In Higgs phase this dual-transformed theory describes a gauge boson of mass et> and a neutral 
Higgs. Gauge coupling e is inversely multiplied to the interaction term between the gauge 
field and the Higgs field, which looks like the strong coupling expansion being done. However, 
when Higgs effects are important, one must take into account the nonpolynomial interaction 
in the Maxwell-like term and the Jacobian in the measure of the Higgs field. Though the 
classical gravity is not affected by the dual transformation, the Jacobian in the measure of 
gravitational field is introduced and this induced Jacobian factor depends on both the gauge 



dynamics and the dimension of spacetime |I3 
The Euler-Lagrange equations are 

jjWwm) = ^wh^ - 11 (3.i3) 

R„u-\g„u(R + 2k) = T^, (3.15) 

where the energy-momentum tensor of dual-transformed theory is 

T° = ^^g p \g^g TK -^;g:)E^E aK + {g^ -\g^^{R p ~d pX ){E a -d aX ) 
+(<?X " \g^9 pC )Q P WM + 9,uV (3.16) 



H — pE^a ds a 5(x a (s a ) 

\/Q 1 J— CO 



X) 



\/9 0= 1 J -°° " s a dS a 

From now on let us look for the counterpart of the vortices in the original theory through 
the derivation of Bogomol'nyi-type bound in the dual-transformed theory. From Eq.(|3.2j) 



and Eq.( [3.7|) , we notice that the classical configurations of both formulations are related by 

\<P\\d,n - eA^) = ^e^d'H". (3.17) 

Together with Eq.( [3.14|) , we obtain 

F* = -gV(H v -d vX )- (3-18) 



For the static counterpart of vortices, Eq. (|3.18|) and the condition that static vortices do not 
carry electric field imply 

F i = H i -K i H o -d iX = 0. (3.19) 

These solitons are characterized by dual electrostatic potential, when N = 1 

$d = fd 2 x^H . (3.20) 

The spatial components of Eq. fl3.14|) are solved with the help of Eq. (|3.19|) , 



K« = ^p§-, (3.21) 



where k d is a constant. Inserting this result into Oi-components of Einstein equations (|3.15|) , 
we obtain 

k d = or ';' -87rGH = C, (3.22) 

where C is an integration constant. Since the vortex solutions of our interest do not satisfy 
the second condition in the sequel, we take the first one from now on. Then T Dl vanishes, 
and the vortices in dual formulation are also spinless as expected. 

The Bogomol'nyi-type bound of the dual-transformed Einstein Maxwell Higgs theory is 
obtained as follows. 



j ' d\^T { 



D 

00 
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d 2 Xy/j 



fir N 2 ~ 



W\ 



1 N 2 ] N ) 

+ -(H 2 + N^HtH^—yidiWd^l + N 2 V + — J2m a 5^ 2 \x - x a ) 

1 Z Vy a=l ) 



d 2 x^H 



2JW 9i { Ho ± \ 



if) fl, F ± -(101- - „') + ^(#o T V2V 



±H,{yw+ e -{\ ( t ) \ 2 -v 2 )) 



(3.23) 



± 



"2^7 



difyiY- 



rJJ 



(■\o\-J v \fl ■ 



2 N 

Z o=l 



To get the last expression we have set iV = 1. Hence, if we choose the scalar potential as 



Eq. (|2.11|) , H becomes 



H = ±-(v 2 -\<P\ 2 ). 



(3.24) 



Since Tq = in this Bogomol'nyi limit and K % i = from Eq. (|3.21|) and Eq. (|3.22|) , ij- 
components of Einstein equations are solved by A = 0. The 00-component of Einstein 
equations is exactly solved as that of the original theory and the solution is Eq.( |2.17|) . 
Substituting this result into Gauss' law (|3.14j) , we obtain a single Bogomol'nyi equation, 
Eq- (CT) . 



IV. Cylindrically Symmetric Solution 

At the outset of our consideration in Einstein Maxwell Higgs model, we had ten second- 
order differential equations (|2.2|)^(|2T4l) (or (|3.13| )^( [3.15| )) for scalar, gauge, and gravitational 
fields even though the system was introduced in (2+1) dimensions after the reduction of a 
dimension along z-axis. However, once we have limited our interest to the static vortex 
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solutions under a specific 4 scalar potential which saturate the Bogomol'nyi-type bound, 
three components of the gauge field and six components of metric tensor have been expressed 
by the scalar field and now there remains only one second-order equation ( 2.19Q to solve. 



In this section let us look for regular and finite-energy vortex solutions of the Bogomol'nyi 
equation of which the base manifold is smooth except for the positions where the point 
particles lie. Let £ be the spatial part of the (2+1) dimensional spacetime. It will be seen 
that there exist solutions such that the space E is a cone or a cylinder asymptotically, or a 
two sphere. 

For the sake of tractability while keeping the main physical properties, let us begin with 
examining the cylindrically symmetric solutions of Bogomol'nyi equation (|2.19|) . The metric 
which is compatible with cylindrically symmetric configurations is of the form 

ds 2 = dt 2 --\-b(r)(dr 2 + r 2 d6 2 ), (4.1) 

lev Y 



where r = ev\fx i x l and < 6 < 2n. Since we consider the regular static solutions of Einstein 
equations, the space £ described by metric b(r) is smooth except for the points where there 
are massive point particles. The global geometry of £ may be characterized by the area 

27j- poo ]_ i-r / J 

A = / drrb(r), the radial distance from the origin p(r) = — / dr Jb(r'), and the 

e 2 v 2 jo ev jo v 

circumference l(r) = — rJbir). 

ev v 

With the aid of gauge transformation, any static cylindrically symmetric field configura- 
tion can be brought into the following ansatz 

= vf(r)e ine = v exp(^ + in6). (4.2) 

Substituting the ansatz into the spatial integral of 00-component of the Einstein equations 
(U) or (ETT5D , we obtain 



itcI d2x ^ 2R 



e 2 jo \\dr ) 4 r 2G(n+M) ^ J 
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where F is the harmonic function factor in Eq. Q2.17 ) which is reduced to a constant in the 



N 



cylindrically symmetric case and M = V] m a is the total mass of point particles which are 

0=1 
now superimposed at the origin. 

For nonsingular solutions, / = \4>\/v should behave as 

/ ~ for n (4.4) 



near the origin. Further conditions are obtained by the requirement that Eq. Q4.3|) is finite. 
Analyzing the behavior of the integrand of Eq. (|4.3|) near the origin by use of Eq. ( |4.4j) , we 
find a condition that 

GM<1, (4.5) 

for finite energy solutions to exist except the trivial one u = 0; analysis in the asymptotic 
region gives the boundary condition for u, 

[ or -00, if < G(n + M) < 1 

u(r = 00) = < _ (4.6) 

I arbitrary number between — 00 and 0, if G(n + M) > 1. 

which is enlarged in compared with that in flat space-time. Now let us examine soliton 
solutions of the Bogomol'nyi equation for the cases of e u ( r=00 )/ 2 ^ and e u ( r =°°)/ 2 — q 
separately. 

(a) e ^ r=0 °)/ 2 ^ 

At first we consider the case that the scalar field e u ' 2 does not vanish at spatial infinity. 
As we see in Eq. fl4.6p , when G(n + M) is smaller than one the boundary value of u at 
r = 00 has to be zero for finite-energy solutions. In this case, it is convenient to introduce a 
variable R such that 

r l-G(n+M) 

R = , ^- . (4.7) 

1 - G{n + M) 
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Then the Bogomol'nyi equation fl2.19|) is rewritten as 



d 2 u dV P ff 1 du 

(4i 



dR? du RdR' 

where V e ff is defined by 

V eff = leM-G{e u -u-l)], (4.9) 

which is an increasing function of u for — oo < u < and has a maximum at u — 0. (See 
Fig. 1.) If we interpret u as a particle position and R as time, Eq. Q4.8| ) is nothing but 
the Newton's equation for a particle of unit mass moving in a potential V e ff and subject to 
a friction. The particle also receives an impact at R = from the delta function term in 
Eq. ( |2.19D . When n = 0, u = is the unique solution and it describes two dimensional flat 



space when there is no particle and a cone when there is a massive particle at the origin. 
When n ^ 0, we now show that there always exists a finite energy solution whose base 
manifold £ is a cone asymptotically. For this we have to show that, if we suitably choose 
the initial parameter Jq i n Eq. ft4.4|) , we can obtain the motion of the hypothetical particle 
such that it starts at negative infinity with the initial velocity given by Eq. ( |4.4|) and stops 
at u = at R = oo. 

First, the behavior of / = e U//2 near the origin is 

fp„G f2G 

f « for n 1 - J JW 2(1 ~ QM) + •••)• (4-10) 



8(1 -CM) 2 

Let R = Ro be an arbitrary large number. If we choose /o sufficiently small, then the higher 
order term in ( |4.10|) can be neglected for R < R Q and the energy of the particle £{Rq) at 
R = Rq is given by 



1 f du 



2 



2 Ui? 



2n 2 



W = obi +y e / /p _ p ^ 7 - *,- . ,•>„,„, • ( 4 - n 



=«o [l _ G(n + M)fRl 

Let us choose R (and / , correspondingly) such that S(R ) < V e ff(0). Then since the 
particle energy decreases as the particle moves due to friction, £ (R) < V e ff(0) for R > Rq. 
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In other words, for sufficiently small /o, the hypothetical particle turns back at a point and 
goes to negative infinity as R —>■ oo. Next, we choose /o arbitrarily large and Rq sufficiently 
small so that (|4.10|) is good at R = i?o- Then it is not difficult to show that for R> i? , 

u(R) > 2lnf r n + ^f(R 2 - R 2 ) + ^R 2 lnR /R, (4.12) 

dV eff 



where VL = max — ^-. The right hand side of Eq. (|4.12|) has the maximum value at 

-oo<u<0 du 

■11/2 



R = R\ 



m + 



In 



^ V' [1-G(n+M)] 



and then u(R\) satisfies, for small R t 



0) 



,„* n /, 4n\l-G(n + M)} \ , , .„ , 

Therefore if we choose /o sufficiently large, then u(Ri) > 0, i.e. the particle goes over 
the hilltop of the potential. From these results, continuity now guarantees the existence 
of the vortex solution connecting the boundary values, u(0) = — oo and u(oo) = 0, for an 
appropriate / . This completes the proof. 

The geometry of S for this solution can be read from the behavior of the metric b(r): it 
is singular at the origin when M^O and the space E has an apex there due to the point 
particle. As r goes to infinity, the radial distance p and the circumference / diverge. The 
asymptotic region in terms of R and 9' (6' = (1 — G(n + M))6) is flat since the solutions 
approach their boundary values exponentially (see Fig. 2-(a),(b)) 

/ w 1 - /ooK ((l - G(n + M))R), (4.14) 

where /oo is a constant determined by the proper behavior of the fields near the origin. 
As shown in Fig. 2-(c), the asymptotic structure of S is a cone with deficit angle 5 = 
2nG(n + M). The solution for M = has been found in Ref. 0. 

Now we consider the case G(n + M) = 1. In this case we define R as R = lnr (— oo < 
R < oo) which reflects the scale symmetry (r — > Xr) of the Bogomol'nyi equation ( |2.19| ) 
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at this critical value, and then it takes the same form as Eq. Q4.8|) with no friction term. 
Therefore in this case the hypothetical particle moves under a conservative force only and 
hence the Bogomol'nyi equation can be integrated to the first order equation 

HSSJ + *•"&=*• (4-15) 

where £ is a constant which is interpreted as the energy of the hypothetical particle which 
is conserved in this case. The particle energy £ is also determined by the initial behavior 
~1 as 

1 / An, \ 2 

= 2n 2 . (4.16) 



£ 2 ^ R ) R- 

In terms of classical mechanics the vortex solution is described as follows: a hypothetical 
particle of unit mass with energy £ = 2n 2 starts at position u = 00 at time R = — 00, climbs 
the hill of the potential V e ff, and finally stops at the top of hill (u = 0) at R = 00. For such 
solutions we must have K//(0) = £ = 2n 2 . From the definition of V e ff, V e ff(0) = F/G and 
the constant F is determined for each n as F = 2Gn 2 . Thus we have completely determined 
free parameters in this case and obtain the solution as the form 

±[- dU = fdR. (4.17) 

\ 1 — e _< ^( eU_u_1 ) 

which is, unfortunately, not integrable to a closed form. It is amusing to note that this kind 
of analysis is not possible in flat case. The behavior of solution near the origin is the same 
as that in Eq.( |4.10|) , so that the space S also has a apex at the origin for M^O solutions. 
However, since b(r) ~ r~ 2 for large r, the radial distance from the origin p diverges but the 



circumference I approaches a finite value I = though r goes to infinity. Then the 

ev 

space S comprises asymptotically a cylinder as shown in Fig. 3. 

Lastly, let us consider the case that G(n + M) is larger than one. If n = 0, GM > 1 
and the only allowed finite energy solution is the trivial one u = as we have seen in Eq. 
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Q4.5|) . Suppose that there exists a solution when G(n + M) > 1 with n^O. According to 
the similar argument given before, the space £ also has an apex if GM is not zero. Since the 
metric b(r) decreases more rapid than 1/r 2 for large r, the radial distance p(r = oo) is finite 
and the circumference l(r = oo) vanishes. Then £ is compact and, since the Euler number 
given in Eq. (|2.14j ) must be nonnegative, two dimensional sphere is the unique candidate. 
From now on let us call the point which corresponds to r = "the south pole" on S 2 and 
that which corresponds to r = oo "the north pole" on S 2 . Then, at the north pole, (f) does 
not vanish and is not well-defined; = 0(r = oo)e me . Therefore there is no regular solution 
in this case. 

(b) e u{r=co)/2 = 

Now let us consider the case that the scalar field e"/ 2 vanishes at r = oo and suppose 
that it behaves like e u l 2 ~ r~ £ for large r. Then the finite energy condition from Eq. (|4.3| ) 
forces e to satisfy G(n + e + M) > 1. Examining the asymptotic behavior of the metric b(r), 
we find that the radial distance p(oo) and the area A of the manifold are finite, and the 
circumference / vanishes at r = oo. Therefore the space S should form a two dimensional 
sphere 5* 2 and it is described by our coordinate (r, 6) except the north pole where a point 
particle may sit; let the mass be M n . The Euler invariant given in Eq. (|2.20|) then should 
be equal to that of S 2 , 

G(n + e + M s + M n )=2, (4.18) 

where M s = M and the subscript s is attached as an indication that it represents the mass 
of particles at the south pole. Now we show that e can not be arbitrary by imposing the 
regularity condition. From the behavior of the metric b(r) at large r, the radial distance 
p{r) behaves near the north pole as p — p(oo) ~ r -G(n+e+M s )+i_ Q n ^e other hand, if we 
choose a coordinate around the north pole rather than around the south pole, we can do 
all the analysis we have done by replacing M s by M n . For example, regularity requires that 
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the scalar field behaves as / ~ s n for s ~ 0, where s is the radial coordinate in the new 
coordinate whose origin is at the north pole; it means that r~ e ~ s n . Also, the radial distance 
near the north pole will behave as p — p(oo) = § Mn+l , Comparing these with Eq. (|4.18|) , 
we get the consistency condition r oc s~ x and thereby e = n. 

Now we discuss the existence of solutions. At first if n = 0, there is no nontrivial solution 
because GM > 1. Next let us consider the case G(n + M s ) ^ 1, i.e. M s ^ M n . Without 
loss of generality we may assume that G(n + M s ) < 1. With r = lni? (— oo < r < oo), the 
Bogomol'nyi equation ([4.8|) is rewritten as 

d2u -e^^IL. (4.19) 



dr 2 du 



Integrating over r from — T to T, 

=T 

2 \dr 



~e 2T V eff \ T T Z T _ T + 2 f dre^V eff . (4.20) 



From the behavior near the poles, / ~ r n (r « 0) and / ~ r £ (r —> oo), it is easy to check 
that the left hand side becomes 

T=T 

■(e 2 -n 2 + 0(e- T )), T -> oo. (4.21) 



1 ( du\ 



2 \dr 



T [1-G(n + M,™ 



On the other hand, the first term of the right hand side in Eq. ( |4.2(J| ) is 0(e~ T ) while the 
second term goes to a positive definite and finite limit as T — ► oo. Therefore, in this case, 
e 2 — n 2 > 0, which is contradictory to the aforementioned condition e — n, i.e. there is no 
regular solution if M s ^ M n . 

The only remaining case is that with M s = M n (= M), or G(n + M) = 1. But in this 
case the second-order Bogomol'nyi equation reduces to the first-order equation (|4.15|) and 



we can simply extend the discussion below Eq. ( (4.15| ). For sphere solutions to exist the 



maximum of the potential V e ff(0) have only to be larger than the energy of the hypothetical 
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particle. [The condition e = n is automatically satisfied.] That is, if F > 2n 2 G for a given 
n, there always exists a unique solution which supports two sphere. If we look at the shape 
of scalar field in r-coordinate, it resembles n/0 nontopological vortex solution: scalar field 
vanishes both at r = and at r = oo. However, since this solution comprises two sphere, 
it can be interpreted as a configuration that two vortices with vorticity n lie both at the 
south pole and at the north pole, and two particles make an apex at each pole if M ^ (see 
Fig. 4). Similar to the cylinder case where F = 2n 2 G, the maximum circumference of two 



2 v / 2rwG 

sphere along the tropical line is which is independent of F. If we regard the point 

ev 



particles as Planck scale strings [T3] parallel to strings of which the symmetry breaking scale 
is lower than the Planck scale, e.g. GUT scale, it may imply a possibility of compactification 
of spatial manifold in the lower symmetry breaking scale. 

V. Conclusion 

In this paper we studied a self-dual system of (2+1) dimensional Einstein Maxwell Higgs 
theory with or without external particles. Bogomol'nyi-type bound for the original and the 
dual-transformed theory has been derived under a specific condition for the form of <p A scalar 
potential. Then, using cylindrical symmetry ansatz, we found all possible soliton solutions of 
the equation. One type of solutions has a Higgs vacuum value as the boundary value and the 
underlying spatial manifold of these solutions is an asymptotic cylinder or a cone. The other, 
which is absent in the flat spacetime case, has a symmetry-restored local maximum value 
as the boundary value and the spatial manifold constitutes two sphere with two vortices or 
vortex-particle composites. These solutions exist when the gravitational constant and the 
mass of external particles satisfy the relation G(n + M) = 1. 

While the existence of cylindrically symmetric solutions of Bogomol'nyi equation has 
been rigorously demonstrated, the stability of those solutions, multi-soliton solutions with 
zero-mode counting, the existence of such solutions away from the Bogomol'nyi limit remain 
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to be clarified. Some aspects of classical self-dual solitons were discussed using the dual- 
transformed form of the model, and the quantum field theoretic issues such as the phase 
transition structure need further study. As the particular choice we made for 4 scalar 
potential was understood by a consideration of supersymmetric models in flat spacetime 
[15J , a supergravity version of the model may be interesting to investigate. 



Note Added 

After writing this paper, we became aware of other references closely related with the 
present paper. The possibility of cosmic strings of cylinder type or 2-sphere type was first 



proposed by Gott jnj. In Abelian Higgs model without point particles, Linet [I7J found 
cylinder and sphere configurations by solving the Bogomol'nyi equations under the choice 
of parameters nG = 1. (But he misinterpreted the sphere solution as representing a vortex- 
antivortex pair; it should be interpreted as representing a vortex-vortex pair as we have seen 
in Sec. 4.) Ortiz |18j studied the same model out of Bogomoln'yi limit including solutions 



with singularity. We thank Ortiz for bringing these references to our attention. 
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Figure Captions 

Figure 1: A typical shape of V e ff(u). V e ff(u) has an asymmetric bell shape with the maxi- 
mum at u — 0. 

Figure 2: Plot of cylindrically symmetric solution with or without point particles at the 
origin for which the underlying space form a cone at asymptotic region with deficit angle 
5 = 2nG(n + M). Parameters chosen in the figures are: M = 0, G = 1/2 and F — 1 (no 
particle); M = 1/2, G = 1/2 and F — 1 (particles at the origin), (a) \<f>\/v vs p, (b) B/b vs 
p and (c) shape of the underlying space E when embedded in three dimensional Euclidean 
space with vertical coordinate denoted as Z. 

Figure 3: Plot of cylindrically symmetric solution with or without point particles at the origin 
for which the underlying space form a cylinder at asymptotic region. Parameters chosen in 
the figures are: M = 0, G = 1 and F = 2 (no particle); M = 1, G = 1/2 and F = 1 
(particles at the origin), (a) \(f>\/v vs p, (b) B/b vs p and (c) shape of the underlying space 
S when embedded in three dimensional Euclidean space with vertical coordinate denoted as 
Z. 

Figure 4: Plot of cylindrically symmetric solution with or without point particles at the 
origin for which the underlying space form a two sphere. Parameters chosen in figures are: 
M — 0, G — 1 and F = 3 (no particle); M — 1, G — 1/2 and F = 3 (particles at the origin). 
(a) \<f>\/v vs p, (b) B/b vs p and (c) shape of the underlying space £ when embedded in three 
dimensional Euclidean space with vertical coordinate denoted as Z. 
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Figure 2-(b) 
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